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Stationary state in N-body System with power law interaction 

Osamu IguchiQ 

Department of Physics, Ochanomizu University, 2-1-1 Ohtuka, Bunkyo, Tokyo,112-8610 Japan 

Many self-gravitating systems often show scaling properties in their mass density, system size, 
velocities and so on. In order to clarify the origin of these scaling properties, we consider the 
stationary state of N-body system with inverse power law interaction. As a simple case, we consider 
the self-similar stationary solution in the collisionless Boltzmann equation with power law potential 
and investigate its stability in terms of a linear symplectic perturbation. The stable scaling solutions 
obtained are characterized by the power index of the potential and the virial ratio of the initial state. 
It is suggested in general that the nonextensive system has much various stable scaling solutions 
than those in the extensive system. 
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I. INTRODUCTION 


There are many self-gravitating systems which are 
characterized by some scaling properties. For example, 
the inter stellar medium shows that its velocity disper¬ 
sion a is power law related with the system size L or 
mass M[jl| (a ~ L 038 ~ M 0 2 ) and isothermal contour 
is characterized by the fractal dimension D ~ 1.36[Q. 
The observations by the Hubble Space Telescope show 
elliptical galaxies have a power law density distribution 
p ~ r~ n (at outer region, n ~ 4 and at inner region, 
n ~ 0.5 — 1.0 for the bright elliptical galaxies and n ~ 2 
for the faint ones||.). The distribution of the galaxies 
and the cluster of galaxies can be characterized by the 
fractal dimension D ~ 2[|J. In cosmological simulations 
based on the standard cold dark matter scenario, the 
density profile shows a power law distribution (at outer 
region, n ~ 3 and at inner region, n ~ 1.0 — 1.5 [|.). 

Recently, in order to study the statistical proper¬ 
ties of a self-gravitating system, we proposed the self- 
gravitating ring model |j, where all particles are mov¬ 
ing, on a circular ring located in three-dimensional space, 
with mutual interaction of gravity in three-dimensional 
space. The numerical simulation shows that the system 
at the intermediate energy scale, where the specific heat 
becomes negative, has some peculiar properties such as 
non-Gaussian and power law velocity distribution^ (n) ~ 
n~ 2 ), scaling mass distribution, and self-similar recurrent 
motion. In this model, the halo particles which belong to 
the intermediate energy scale are considered to play an 
important role in realizing such specific characters. 


We are interested in the origin of these scaling prop¬ 
erties from statistical mechanical point of view. In or¬ 
der to study the statistical properties of long range in¬ 
teraction such as gravity, Ising model, and spin glass, 
the model with power law potential has been used and 
revealed anomalous properties || |], 10 . For exam¬ 
ple, a gravitational-like phase transitionH, reduction of 
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mixin g|, and long relaxation |I(J are observed. Using 
a model with an attractive 1 jr a potential in general D- 
dimensional space, we can control the extensivity of the 
system and the specific heat by changing the spatial di¬ 
mension D and the exponent of inverse power of the po¬ 
tential a. 

In this paper, we study the quasi-equilibrium state of 
N-body system with a power law potential. As a first 
step, we consider the collisionless Boltzmann equation 
(CBE) in place of N-body system and derive the self¬ 
similar stationary solution of CBE which has a scaling 
property appearing in the quasi-equilibrium state and 
discuss the linear stability by the use of energy functional 
analysis Jll|, 0, |I|, |l §]. 

In section [H], we show some general properties of N- 
body systems with power law potential. In section III, 
we derive the self-similar stationary solution of CBE with 
an attractive l/r“ potential assuming spherical symme¬ 
try and isotropic orbit case in D-dimensional space. Sta¬ 
bility for the linear perturbation around the self-similar 
stationary solution is investigated in section IV. Section 
^ is devoted to discussion. 


II. N-BODY SYSTEM WITH POWER LAW 
POTENTIAL 


In this section, we show some general characteristic 
properties of N-body system with power law potential. 

We write the Hamiltonian for the N-body system with 
power law potential in the form: 


JV 


^ = Gm 2 


i=1 


2 m 


i<j 


(1) 


where rij := r 7 ; — r 7 - and a controls the range of inter¬ 
action. 

In this system, the virial equilibrium condition be¬ 
comes 


2 < K > +a < $ >= 0, (2) 

where < K > is the time averaged kinetic energy and 
< $ > is the time averaged potential energy. From the 
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TABLE I: In the system with an attractive 1 jr a potential in 
D-dimensional space, the property of the specific heat and 
the extensivity is shown. 



0 < a < 2 

a > 2 

specific heat 

negative 

positive 


a < D 

a > D 

extensivity 

nonextensive extensive 


expression H = K + we have 



a 


< K >= 


2 -a 
2 


< $ > . 


( 3 ) 


From Eq.(||), the signature of the specific heat is deter¬ 
mined by the signature of the term —(2 — a)/a. 

In order to clarify the extensivity / nonextensivity of 
the system, we focus on the N dependence of the poten¬ 
tial energy $ under fixing the number density Al/L D ]l5[: 

N a/n) 

drr 0 - 1 ^ 01 ~ N l ~ a/D . (4) 



If the N dependence of the potential energy per particle 
disappears for N —» cxd. we define the system extensive. 
Otherwise, we define the system nonextensive. In the 
case of the gravity in D-dimensional space, since a = 
D — 2, the system is always nonextensive. We summarize 
the signature of the specific heat of the system and the 
extensivity in Table .[j| 


III. SELF-SIMILAR STATIONARY SOLUTION 
IN COLLISIONLESS BOLTZMANN EQUATION 
(CBE) 


In this section, we derive a self-similar stationary so¬ 
lution in the collisionless Boltzmann equation (CBE); 

< 5 > 

where / = f{x,v,t) is a mass distribution function and 
[. A, B] denotes the Poisson bracket. 

The stationary solution f 0 satisfies the following equa¬ 
tion, 




^ dt dw l 


( 6 ) 


where w l = {x, d}. 

For the coupled CBE and Poisson equation, R.N. Hen- 
riksen and L.M. Widrow|l6|] studied the self-similar 
stationary solution in CBE with spherical symmetry 
in three-dimensional space by applying the systematic 
method which is based on the work of B. Carter and 
R.N. Henriksen[[l7j. 


Following R.N. Henriksen and L.M. Widrowfl(|. we 
study the self-similar stationary solution with spheri¬ 
cal symmetry and isotropic orbit case in general D- 
dimensional space. The case that D = 3 and 

a = 1 reduces to the work by R.N. Henriksen and 
L.M. Widrow[|Hi|. By the generalization of the spatial 
dimension D and the exponent of power of potential a, 
we intend to investigate the relation between the exten¬ 
sivity of the system and the self-similarity. 

From Eq.(|(j)), mass distribution function f(r,v) obeys 

vd r f - d r $d v f = 0, (7) 

where v := + Vg + and is a potential. The 

potential $ satisfies the following equation, 

dr (r a+1 9 r 4>) = S 2 d G J v D ~ x fdv , (8) 

where Sd ■= 2tt d ^ 2 /T(D/2). In the case of a = D — 2, 
the above equation corresponds to Poisson equation. 

A self-similar stationary solution satisfies the following 
equation, 


Ckf = 0, (9) 

where 

Ck '■= k z di = 6rd r + wd v + nmd m ( 10 ) 

is a Lie derivative with respect to the vector k in phase 
space, and 5 , u, and /i are arbitrary constants. 

In a dimensional space of length, velocity, and mass, 
we introduce vectors a = (<5, is, /i) and df. The vector 
a = (<5, u, fx) describes changes in the logarithms of di¬ 
mensional quantities. Each dimensional quantity / in 
the problem has its dimension represented by the vector 
df. Using these vectors a and df, the action of k reads 

Ckf = {d f ■ a)f. (11) 

The dimensional quantities in current problem /, <F, 
and G have the following dimensional covectors, 

d f = {—D, —D, 1), 

= ( 0 , 2 , 0 ), ( 12 ) 

da = («, 2 ,- 1 ). 

The requirement of the invariance of G under the rescal¬ 
ing group action (fio]) implies da ■ a = 0, 

/x = ad + 2v. (13) 

The dimensional space is reduced to the subspace of 
(length, velocity), wherein the rescaling group element 
a = ( 5 , ix) and 


df = (a — D,2— D), 

= ( 0 , 2 ). 


( 14 ) 
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Here we define the new coordinate R(r) and X in re¬ 
placement of the original coordinate r and v such that 


CkR = 1 , 
C k X = 0. 


From Eqs.(|l5]) and ([il]), we choose 


JR 


i<5| = 
v = Xe 1 


iyR 


(15) 

(16) 


(17) 

(18) 


The transformation from the original coordinate (r, v) 
to the self-similar coordinate (R,X) is shown in Ap¬ 
pendix]^. 

Under the new coordinate, these physical quantities / 
and $ can be written in the form: 

f(X,R) = J(x) e - [iD - a)s+{D - 2 W R , (19) 

$(X,R) = ${X)e 2vR . (20) 


Substituting Eqs.(|l9|) and (20) into Eqs.(^) and 
these equations for a bounded solution yield 


din/ _ [D — 2 + (D — a)^\ X 2 


d In X 

2 v 2 \ 8\ D ~ a ~ 2 


X 2 
u5 


2 $ 


$ = S 2 D G 


- 2 $ 
X D ~ l 


( 21 ) 

IdX. (22) 


Without loss of generality, we can set v = 1. 

Solving Eqs.([ll|) and (|22|), we have the following solu¬ 
tion, 


f = C\X 2 + 2<f>|~[(' D_Q )' 5+ ( I3_2 )]/ 2 , 


(23) 


where 


C = 


, 5 \ |< 5 | 


D-a-2 


T{D/2)Y{2 + (a - D)S/2) 


2n D G\ - 2$|(“- £l )U2r([4 — D + (a — D)S — D}/2) 


(24) 


and the following condition must be satisfied, 

(D - a)S < 4 - D. (25) 

Since 4> < 0, from Eq. (p2|) we obtain the additional con¬ 
dition, 


a8 < -2. (26) 

If these condition Eqs.(p5|) and (^) are actually sat¬ 
isfied, we have the bounded self-similar stationary solu¬ 
tion ( p3| ) and (p4|). The mass distribution function /, the 
mass density p, and the velocity distribution f(v) become 
respectively 

f(r, v) = c\2E\- [{D - a)6+{D - 2 ^/ 2 , (27) 

p := S D J dw D - 1 /(r,u)~r a - D+2 / i , (28) 

f(v) := S D J drr D - 1 f(r,v)~v a5+2 ~ D , (29) 


where E denotes the mean field energy: 

E := i v 2 + $o- (30) 

The ratio of the average of the kinetic energy to the 
potential energy is as follows: 


< 4>o > _ (D - a)6 - 4 

< K > = 2D 


Since the solution (^) we obtained is a bounded solution, 
which satisfies the following condition, 


(.D — a)S — 4 
2D 


< - 1 , 


(32) 


the specific heat of the self-similar stationary solution is 
always negative. The d* corresponds to a virial equilib¬ 
rium state: 


J-i (* 7 ^) 

1 arbitrary (a = D). 


(33) 


If (D — a)(5 — 5*) < 0, the potential energy in this state 
is more dominant than in the virialized state. 

The relation between pressure P and mass density p 
can be written in the form: 


P 


i+i+^d 


D)5/2 


(34) 


The above equation of state corresponds to Polytropes 
gas when identifying the Polytropes index n as 1 + (a — 
D)S/2. Note that for a = D, there is no self-similar 
stationary solution corresponding to the isothermal state. 

As for gravity case (a = D — 2), the above solution ( |24| ) 
and (|7|) in D = 3 corresponds to the solution derived by 
R.N. Henriksen and L.M. Widrow|l(ij|. For D = 1 and 
D = 2 where a = D — 2 < 0, we show the self-similar 
stationary solution in Appendix 


IV. LINEAR PERTURBATION ANALYSIS 

In this section, we investigate the stability of the self¬ 
similar stationary solution derived in the previous section 
for a symplectic linear perturbation, by energy functional 
analysis Ju]. @- 

As for the linear stability of the stationary solution 
in CBE of the gravity in three-dimensional space, there 
has been much research 0 0 0 0 0 0- 0- 0, 
|22| , |23| , p4| . For the stationary state, assuming spherical 
symmetry, characterized by the mass distribution func¬ 
tion /o specified as a function of the mean held energy E 
and the squared angular momentum J 2 , if dfo/dE < 0 
and dfo/dJ 2 < 0, then the system is stable to the linear 
perturbation. 

Following the work by J. Perez and J.J. Aly[^3| where 
the stability of stationary solution in the coupled CBE 
and Poisson equation with spherical symmetry in three- 
dimensional space was studied, we study the stability of 
the solution obtained in the previous section. 
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First, we explain a symplectic linear perturbation by 
energy functional analysis ifTH |ij| pj], |ix|. In term of the 
mass distribution function f(x,v,t), the Hamiltonian H 
is written as follows, 

H = J dT^-f{x,v,t) 

+ ^ J dT JdT'Q(\x — x'\)f(x, v, v', t),(35) 

where dr := d D xd D v is the 2D-dimensional phase vol¬ 
ume element and the kernel Q satisfies 

$(*) = G J dT'G(\x - x'\)f(x',v',t). (36) 

We consider a small perturbation around the station¬ 
ary solution / 0 . The distribution function and Hamilto¬ 
nian can be expanded around the stationary solution as 
follows, 

= fo + <5 (1) / + d (2) / H-, (37) 

H = H 0 + 6 {1) H + 6 {2) H + ■■■ . (38) 

Here we consider any symplectic perturbation, which 
can be generated from the stationary solution fo by use 
of a canonical transformation. By using some generating 
function K , any symplectic deformation can be expressed 
in the form: 

/ = e [if ’ ] /o- (39) 

From the above definition (|39|), / can also be expressed 
as follows, 

/= fo + [K,f 0 ] + ^[K,[K,fo]} 

+l i [K,[K,[KJ 0 ]] + ... . ( 40 ) 

Introducing a small parameter e which represents the 
amplitude of the perturbation, we expand the generating 
function K as 


K = eK w + e 2 A' (2) + e 3 if (3) + • • • . (41) 

Further identifying = e n K^ n \ we obtain the per¬ 
turbed quantities in the Eq. (|37|) in the form: 


s {1) f = [g (1) , fo], 

(42) 

t-h | CN 

+ 

c? 

II 

CM 

(43) 

The first order term in Eq.(l3^) becomes 


S^H = J dTE[g^,f 0 }, 

(44) 

where E is the energy of a particle, 


V 2 

E:=— + <f>o, 

(45) 


where <f>o is the potential energy generated by fo- Since 
E and fo are conserved quantities, 5^H = 0. 

The next order term in Eq. (|38|) yields 

S (2) H = J dTE[g^\fo} + \ J drE[g^,[g^Jo]} 

+ jJ dT J dT'Q(\x ~ x'\)[g {1 \ fo][g (1) ', fo\. (46) 

The first term in Eq.([l6|) also vanishes and by the inte¬ 
gration by parts, Eq. (|4q) is rewritten in the form: 

6 {2)h = -\ J dT[g^,f 0 ){g^,E] 

+ f J dT j rfT'6?(|* - * # |)[ff Cl) , /o][ff (1)/ , /a- (47) 

Hereafter we consider the case that the stationary so¬ 
lution fo is a function of only the energy E. In this case, 
we obtain 

[3 (1) ,/o] = F E [gW,E\, (48) 

Jd D v[gU,f 0 ] = J d D vd x (F E vg^), (49) 

where F E := d E fo- 

Integrating by parts and using Eqs.([l8|) and (|4g|), we 
have 

S^H = 

\ J dT(-F E )\lgW,E]\ 2 + ^ J dTd x {-F E vg «) 

x J dY'd x ’{—F' E v'g^ lS> ')Q{\x — x'\). (50) 

The linear perturbation g W has two kinds of gauge 
mode, (a) <?W = g^\E). In this case, the linear pertur¬ 
bation of the mass distribution f is trivially zero. ( b) 
gT) = av (a is a constant.). This perturbation means 
the translation of the center of mass. In order to con¬ 
sider the physical perturbation, we investigate the linear 
perturbation excluding the above gauge modes. 

The stability for the linear perturbation (2^, |2l| reads: 

If 6&H > 0, then the system is stable. (51) 

A. spherical mode 

Since the first order perturbed potential (5^)$(r) sat¬ 
isfies 

^d r (r a+1 d r 5^<S>(r)) =S D G J d D v'[g^',Q 

= SnG-^dr (V 0 - 1 J d D v'F’ E v r 'g , (52) 

the spatial derivative of j( 1 )$(r) becomes 

J d D v'F’ E v r 'g < 


(53) 
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From Eqs.([50|) and (|53|), we have 
2 5 {2) H 

= jdT(-F E )\\g (1) , E}\ 2 
-j dTd r (-F E v r g (1) )5 < ' 1) <f>(r) 

= J dT(-F E )\[g^,E]\ 2 -S D G J 

x J d D v(—F E vg^) j d D v'(—F E v r g ^ ). (54) 

Introducing new variables, 

9 {1) =: rv r p(r,v,t), (55) 

and using Schwartz’s inequality, we have 


nonextensive extensive 



FIG. 1: Stability chart for the one-dimensional case (D = 1). 
The dark region corresponds to the stable region (163) in the 
parameter space (5, a). The property of the specific heat and 
the extensivity also are shown. 


2 S^H = J dT(-F E )\[prv r ,E}\ 2 -GS D J ■ 


d D x 

,a-D +2 


x J d D v[-F E r(v r ) 2 p] J d D v'[-F’ E r{v r ) 2 p’\ 
> J dT(-F E )\[prv r ,E]\ 2 -GS D J 

x J d D v[-F E r(v r ) V] j d D v'[-F' E r{v r ) 2 } 
= Jdn-F E )[\[grF,E]\ 2 ^ GS ^ Vr)2 ^^ 


r a-D +2 f ’ 


where po is non-perturbed mass density: 

Po := j d D vf 0 = j d D v{-F E ){v r f. (57) 

Using the property of the Poisson bracket, and the fact 
that the integral of the Poisson bracket over the phase 
space vanishes, the equation ( |56| ) can be rewritten in the 
form: 


2 S^H > J dT(-F E )^\[prv r ,E}\ 2 - 


2 GS E {rv r ) 2 p 2 po 


y*ol — D-\-2 


= J dT(-F E ){(rv r ) 2 \[p,E}\ 2 + \p\ 2 \[rv r ,E}\ 2 

GS D (rv r ) 2 p 2 p Q 


+rv r [p 1 E)[rv r ,E] — 


y%ol — D-\-2 


= J dT(-F E ){(rv r ) 2 \[p,E]\ 2 + \p\ 2 \[rv r ,E]\ 2 
2 rv r [rv r ,E], 

-\p\*\[rv r ,E]\ 2 - 


+ [p z rv r [rv r ,E],E} - \p\ 2 rv r [[rv r , E], E] 
GS D (rv r ) 2 p 2 po 


r>QL -D + 2 


= j dT(-F E ){(rv r ) 2 \[p,E}\ 2 


i 12 r rr r pi pi GS D (rv r ) 2 p 2 po , , . 

-H rV ll rV , E }’ E \ - -a-D+2 - f '( 58 ) 


Using the following relation 
/d 2 <h 

[rv r ,E\,E\ = —rv r 


o , 3d$o 
r dr 


dr 2 

f GSdPo 2 - a #q 


1 rpCX. — D -(-2 


dr 


,(59) 


we obtain the final expression in the form: 

5 {2) H > 

|2 2 - a o 
dr 


(56) \ J dV{-F E ){rv r ) 2 [\{p,E}\ 2 + \pf 


(60) 

From Eqs.@), (|l7|), and (|60|), if S^H = 0 when F E < 0 
and a < 2, 1 = 0. Since this is a gauge mode, we 
conclude that 

If F e < 0 and a < 2, then S^H > 0. (61) 

From the self-similar stationary solution Eq.(^3j), we 
have 

F e = sgn (E)[(a - D)S + 2 - D]C\2E\^ a - D "> s - D ^ 2 .(Q2) 

As an explicit example, we consider D = 1 case. From 
Eqs.(|||), (pti|), and (|6l|), if the following condition 
is satisfied, the self-similar stationary solution Eq.(27) is 
stable. 


1 r 2 

S<-~ 


(1 < a < 2) 

(0 < a < 1) (63) 


In Fig.|l[ we show the region where exits the stable self¬ 
similar stationary solution in the parameter space (a, d). 

Note that in the above calculation, we use the integra¬ 
tion by parts and neglect the surface term at the origin. 
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Since the self-similar stationary solution obtained in this 
paper is singular at the boundary, the surface term can 
not be neglected in general. In the realistic situation, 
however, the self-similarity appears in the intermediate 
scale since the system has a cut off scale in the short 
distance. We suppose that the self-similar stationary so¬ 
lution can be connected with some regular solution near 
the boundary by regularization such as r —> ( r 2 + a 2 ) 1 / 2 
where a is a cut off scale, and the boundary term can be 
neglected. 


we can rewrite the equation (68) in the form: 

6&H 


dr 


-Fe 


+ ys^> j dDx y r$ °) 2 J d D v{-F E )\w 

-\w\ 2 d r <r a w 2 d r <r 0 >. (70) 


B. aspherical mode 


Next, we consider the aspherical mode. Since it is dif¬ 
ficult to analyze a general case, we study the gravity case 
in D-dimensional space (a = D — 2). 

By the integral of Poisson equation over the configu¬ 
ration space and integration by parts, we have 


By straightforward calculation, we obtain 


V 2 a r $ 0 = S D J d D vF E d r % + (D — 1)^. (71) 

By using Wirtinger’s inequality, we have 



D - 1 
^2 



dVL > 0, 


(72) 


1 


where 


J d D x\VS^^\ 2 = J dTSWfSWQ, (64) 


:= j d D vQ8 (1) f. 


(65) 

From Eqs.(ff7P, (|64|), and (j35j), we have 

S^H =\J dr [ FFF - J- | d D x\VS^^\ 2 . (66) 

Here we introduce 5^f as follows, 

5 {1) f =: F E 6 w $ + S w f. (67) 


where V s := V — r-r0-. 

° _ | r | o r 

From Eqs. ©, (ffll), and (^|]), we get the final expres¬ 
sion in the form: 

5™H 


+ 


2 J -F e 
1 


2 S D 


d u x(d r ^oY 


\d r 


|V S 


D- 1, 


i r ntfWf ) 2 i 
~ 2 J -F e 2 S D 


J d D x(d r <ry\d r w\ 2 . (73) 


From Eq.(|7^), if S^H = 0 when Fe < 0, then 5^ f = 
0 or </'> = av. Since this is a gauge mode, we conclude 
that 


If F e < 0, then H > 0. 


(74) 


Substituting Eq.(|67|) into Eq. (|66|) , we obtain 

S^H 


dr 


+ f e { ( 5 (1) $) 2 - 2<5 (1) /<5 (1) 


$ 


G 

~2SY>. 

-S D 


dr 


d D cc|V<5 (1) $| : 

(<5«/) 2 


—F e 2 Sd 
d D v(-F E ) |<^$| 2 


J d D a;||V5 (1) $| 


( 68 ) 


Moreover, using the new variable w which is defined 

by 


=: w(x, t)d r $o, (69) 


This condition is weaker than the condition ( pl|). I n the 
gravity case (a = D — 2), from Eqs.(|25|), (|26|), (J32|), and 
& if the following condition is satisfied, the self-similar 
stationary solution Eq.(^) is stable. 


<5* 


<5 <-1^2 (2<D<2 + V2) 

5 <2 — D (2 + y/2 < D < A). K ’ 


In Fig.|^, we show the region where exists the stable self¬ 
similar stationary solution in the parameter space ( D , <5). 
This stability condition ( pl| ) is consistent with the work 
by J. Perez and J.J. Alyp3[(q = 1, D = 3). 


V. DISCUSSION 

We studied the self-similar stationary solution in colli¬ 
sionless Boltzmann equation with the attractive l/r a po¬ 
tential. Assuming the spherical symmetric and isotropic 
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nonextensive 



teraction such as gravity. In the realistic situation, since 
the anisotropic velocity space is important, we would like 
to extend this analysis to the anisotropic case in our fu¬ 
ture work. 
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FIG. 2: Stability chart for the gravity case (a = D — 2). 
The dark region corresponds to the stable region (|75| ) in the 
parameter space ($, D). The property of the specific heat and 
the extensivity also are shown. 


orbit in D-dinrensional space, we investigate the linear 
stability of the solution. In the above model, we can 
control the extensivity of the system and the signature 
of the specific heat by changing the spatial dimension D 
and the exponent of inverse power of the potential a. 

The self-similar stationary solution can be expressed in 
the form of the power law of the energy. The exponent of 
the power is determined by the power of the potential a, 
spatial dimension U, and the scaling parameter 5. Here 
we interpret S as a parameter which denotes the virial 
ratio of the initial state. 

By use of the energy functional approach, we investi¬ 
gated the stability of the self-similar stationary solution 
in terms of a symplectic linear perturbation. As for the 
spherical symmetric and isotropic orbit of the gravity in 
D-dimensional space (a = D — 2), we found that the sys¬ 
tem is stable if the mass distribution function decreases 
monotonically and the spatial dimension is less than 4. 
As for the power-law potential in one-dimensional space 
(D = 1), we found that the system is stable if the mass 
distribution function decreases monotonically and the in¬ 
verse power index of thepotential is less than 2. The self- 
gravitating ring model Q is similar to the case of a = 1 
in one-dimensional space. From the form of the velocity 
distribution obtained by a numerical simulation, <5 ~ —3. 
This case belongs to the stable self-similar stationary so¬ 
lution. 

The stable self-similar stationary solution we obtained 
includes the virial equilibrium state in the case of 5 = <5* 
As for the extensivity of the system, the nonextensive 
system has far more stable scaling solutions than the ex¬ 
tensive system in the parameter space (<5, a, D). 

In the time evolution of the collisionless system, assum¬ 
ing the spherical symmetry and isothermal case, Larson- 
Penston solution which shows self-similar collapse is the 
attractor[^6|. By such a self-similar time evolution of 
system, we expect that the class of the stable self-similar 
stationary solution obtained in this paper plays an impor¬ 
tant role as a quasi-equilibrium state with a long range in- 


APPENDIX A: TRANSFORMATION TO 
SELF-SIMILAR COORDINATE 

Using the original coordinate (r, v), the self-similar co¬ 
ordinate (R, X) is defined by 

R := <r 1 ln(r|<y|), (Al) 

X := ve~ vR . (A2) 

The derivative with respect to the original coordinate 
can be expressed by the self-similar coordinate in the 
form: 


d r 

d v 


dR 

dr 
sgn((5)e 
dX 
dv 
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x 

-SR 


dX 

dr 


dx, 


(Or - vXdx) 


dx = 


-vR 


dx- 


(A3) 

(A4) 


APPENDIX B: STABILITY CONDITION FOR 
GRAVITY CASE IN D = 1 AND 2 (a = D - 2) 


For the case that the potential •!> is positive, Eq.(|2^) 
is modified as 


2v z 


2 + (£>- 2 )- 


$ = S 2 d G X u ~ l fd,X. (Bl) 


1. D = 1 case 


Since the potential $ is positive, from Eq.(Bl), 

6 < 2. (B2) 


By integrating Eq.([Bl|), we have a self-similar solution 
(E 3 I) and 


C = 


| 2 -( 5 | r (< 5 - 1 / 2 )| 2$|' 5 


V^GT(S - 1 ) 

if the following condition is satisfied: 

5 > 1. 


(B3) 


(B4) 
















From Eqs. ( |B2| ), (|B4|), ( |6l| ) and (|7J), the stability 
condition for linear perturbation yields 


and the integral constant of ([230 is 


1 < 5 < 2. (B5) 

The ratio of the average of the kinetic energy to the 
potential energy is the same as Eq.(|3l|) in D = 1. How¬ 
ever, if S < 2, the integral of the kinetic energy over the 
velocity space diverges. For this reason, there do not 
exist a stable self-similar stationary solution in D = 1. 


2r(<5)|2<J>|' 5 
7t 5 / 2 GT(<5 — 1/2) ’ 


(B7) 


From Eqs. (61), (|74|), and (B6), the stability condition 
against linear perturbation yields 


2. D = 2 case 


If the potential $ is negative, the condition that 
the bounded self-similar solution exists is the same as 
Eqs.(^5|) and (^6|). Since this case does not satisfy the 
condition (^gJ) , the only case is that $ > 0. 

In this case, from Eq.(Bl), the condition that a 
bounded self-similar solution exists yields 


6 > 

2 


(B6) 


5 > 


1 

2 ' 


(B8) 


The ratio of the average of the kinetic energy to the po¬ 
tential energy is the same as Eq. © in D = 2. However, 
similar to the D = 1 case, if 6 < 2, the integral of the 
kinetic energy over the velocity space diverges. Finally, 
if S > 2, the self-similar stationary solution in D = 2 is 
stable. In this case, the specific heat is always positive. 
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